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We have studied the interplay between disorder and cooperative scattering for single scat- 
tering limit in the presence of a driving laser. Analytical results have been derived and we 
have observed cooperative scattering effects in a variety of experiments, ranging from ther- 
mal atoms in an optical dipole trap, atoms released from a dark MOT and atoms in a BEG, 
consistent with our theoretical predictions. 
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1. Introduction 

The interaction of quasi resonant light with a large cloud of atoms has been studied 
for many years, starting with seminal work by Dicke ^ and with renewed interest 
in the context of entanglement which such systems are expected to contain. Often 
continuous density distributions have been assumed which allow for analytical ex- 
pressions to be obtained. The role of fluctuations of the atomic density is however 
by itself at the origin of interesting phenomena, such as Anderson localisation of 
light (0 [2) ■ In a series of theoretical and experimental studies, we have recently 
addressed the question of the quasi-resonant interaction of light with clouds of cold 
atoms, bridging the gap from single atom behavior, effects dominated by disorder 
to a mean field regime, where a continuous density distribution is the relevant 
description. 

In this paper, we present a theoretical model we use to describe the collective 
atomic response under continuous excitation of a low intensity light field, with 
additional details compared to Q- We then present experiments which have been 
performed in complement to those reported in 0, using atoms in a dipole trap 
as well as atoms in a magnetic trap, both above and below the Bose-Einstein 
condensation temperature. 
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2. Theoretical Description 



Hamiltonian and state of the system ~ We consider a cloud of N two level atoms 
(positions r^, lower and upper states \gj) and \ej) respectively, transition frequency 
Wq, excited state lifetime l/F), excited by a quasi-resonant incident laser propagat- 
ing along the direction (wave vector ko) and with frequency wq = f^a + Aq. The 
atom-field interaction Hamiltonian is in the rotating-wave approximation (RWA) 

^ ro 1 ^ 

H = hY, ^a.e*^"*-^'^-'-' + h.c. +;i^^[5fca,4e'^'=*-'''"^ +h.c.] . (1) 



j=i k 



Here, is the Rabi frequency of the interaction between an atom and the classical 
pump mode, dj = \gj){ej\ is the lowering operator for atom j, Ok is the photon 
annihilation operator, and = d^Ju:k/ {heoVph) describes the coupling between 
the atom and the vacuum modes with volume Vph and frequency cj^ = + A^. 
We assume that all atoms are driven by the unperturbed incident laser beam, thus 
neglecting dephasing by atoms along the laser path or by near field effects, which 
could arise for large spatial densities. Calling |0)a = |5i, ■■,gN) the atomic ground 
state and |j)a = \9i, ej, Bn) the state where only the atom j is excited, we 
assume that the total state of the system has the following form (0, : 

N 

|^(t)) = a(t)|0).|0)k + e-^«* J]/3,(t)e^''«-^ |j%|0)k + j;7k(t)|0)a|l)k • (2) 

j=i k 

The above expression assumes that only states with at most one atomic excitation 
contribute to the effects here described (fiOf) . 

Time evolution of the system - The time evolution of the amplitudes is obtained 
by inserting the Hamiltonian ([T]) and the ansatz ^ into the Schrodinger equation, 
dtl^it)) = -{i/h)H\^{t)): 

N 

d = -ifio^/3, , (3) 
4- = ^ -i^^a-iY, 5fc7ke^(^«-'''=)*+*('^-''")"^ , (4) 



7k = - igke- 



N 

-i(Ao-Ak)tSr^ 5 -i(k-ko)-r, 



^Pje-'^^-^''^'^. (5) 



Integrating Eq. ^ over time and substituting 7k (i) in Eq. Q we obtain: 

N ,.t 
i_ 1 



k m=l 



Assuming the Markov approximation (valid for r^r o"r/c, where cJr is the size of 
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the atomic cloud and rjv is the cooperative decay time), we can approximate 



Then, going to continuous momentum space via — )■ Vph{2'iT)^^ dkk"^ J dfik 
(where dfl^ = sin 9 d9 d(j)) and neglecting saturation assuming a ~ 1, we obtain 

N 

= iAo^j -^^O-l^Y^ JjmPm , (8) 
m=l 

where T = {Vph/TTc)k'^gl^ and (jH]) 

^ =J_ [ rfj^^g^(k-ko).(r,-r„) ^ ^-.ko.(.,-.,„) sm(fco|r,-r^|) 
47r 7 ko\rj - r^l 

Due to the presence of the driving term, the solution will evolve quickly toward 
the driven timed Dicke state Q, characterized by 

m = (10) 



Once inserted the ansatz (10), Eq. ([8|) yields 

P = -^VNno+{iAo-^TNsN)P , (11) 

where 

s^ = — d(l) desme\SN{ko,9,(l))\'^ . (12) 
47r Jo Jo 

and S'Ar(k) = (1/iV) XljLi exp[— i(k — kg) • fj] is the structure factor of the atomic 
cloud. In steady state we find 



« (13) 

In fig jTI w e compare the probability to find atoms in the driven timed Dicked 
state (10), PoTDit) = |/3(i)P (red continuous line) obtained from Eq. (11), and 
the probability that atoms are excited, Pe{t) = XljLi (blue dashed line) 

obtained solving numerically Eq.([8]), for a spherical gaussian cloud with N = 4000 
atoms, size a = koCTr = 10 and a pump beam with Aq = lOE and VLq = O.IF. We 
observe that the exact state tends toward the driven timed Dicke state. 

Forces in the Markov approximation. — The two terms in the Hamiltonian ([T]) 
yield two different contributions to radiation pressure force: 

Fa, +Fei = -Vr,i^ . (14) 

We will be interested in the average absorption force, = '^ji^aj)^ and emis- 
sion force, Fe = jj^j{^ej) ^ acting on the center of mass of the whole cloud, 
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Figure 1. (color online) PoTD^t) = |/J(t)P (red solid line) and = |/3j(t)p (blue dashed line) as a 



function of Ft obtained using the evolution equations (ITTf and (jsl, respectively. The simulation parameters 
are N = 4000, cr = 10 for a gaussian spherical cloud, Aq = lOF ai 



:tively. 
"and Qo = 0.1 F. 



Ffj + Fe = macM, where acM is the center-of-mass acceleration and m the mass of 
one atom. The first term, Faj = |^kof^o['5"je*^°*~*''°'''' — h.c] results from the recoil 
received upon absorption of a photon from the pump laser and has an expectation 



value on the timed Dicke state ( 10 ) given by: 



Im[/3(i)] 



(15) 



1. The second contribution, Fgj 



where we assumed again a " 

i 7ik(7fc[(Tjd|^e*^''*~*'^''"^ — h.c], results from the emission of a photon into any 
direction k. The expectation value on the general state ^ is: 



(16) 



Substituting the time integral of 7k (i) from Eq. ([s]) and inserting the timed Dicke 
state from Eq. ( 10 ) we obtain the average emission force: 



k 



Pit) / dtV('^^-^«)*'/3*(t-t') + c.c. 

^0 



(17) 



In the Markov approximation and going to continuous momentum space we find 



(18) 



where 



In 



2n 



d0 sin6'cos0 \SNiko, 0, < 



(19) 



Finally, using Eq. (13) in Eqs. (15) and (18), the average steady-state radiation 



force acting on the center of mass of the atomic cloud is 



F„ + Fe = ^koT 



4A2 + N^r^sl 



(sN - In)- 



(20) 
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The common prefactor can be obtained from the standard low saturation single- 
atom radiation force Fi = KkoT^Q/ (4AQ+r^) by substituting the natural linewidth 
by the collective linewidth, F — )• NTs^, and the Rabi frequency by the collective 
Rabi frequency, Qq — ^ VN^Iq. Additionally, the cooperative radiation pressure 
force is weighted by the difference of structure factors, sn — fw, where the sn 
part corresponds to the cooperative absorption process and the fjy part to the 
cooperative emission. For smooth density distributions n(r), one could compute 
the structure functions by replacing the sum with an integral (sat — )• Soo and 
/tv — )• foo)- However, we have shown (jj that in this way we miss the role of the 
disorder in the scattering process and the crossover from single atom radiation 
force and cooperative scattering. Instead, estimating the fluctuations of sat and /at 
one finds that (see figj2t 



1 



In ~ /o 



(21) 



atom force is 



Using (21), the ratio between the cooperative radiation force (20) and the single 

[1 + N{soo - foo)] (22) 



Fr. 



4A2 + F2 



Fi 4A2 + F2(l + Ars^)2 



Assuming a smooth Gaussian density distribution with ellipsoidal shape, 
noexp[— (x^ + y'^)/2a'l — z^/2a1], the structure factor is 5oo(A;o, ^, </>) = 
exp{— cr^[sin2 + 'rf{cos9 — l)^]/2}, where a = k^ar and r] = azj^r is the aspect 
ratio. For elongated clouds, ?/ > 1, 



fiv) 

I oo 



erf 



f7(2r/2 



erf 



„2 (r,) 



4(t2 



(23) 



For spherical clouds (rj = 1) and for a ^> 1 one finds Sqo ~ l/4a"2 and Soo — fa 



1/8(T^. For cj, 7? S> 1, can be approximated by s'^' ~ Sooy/T^Fe^^ [1 — erf(F)], 
where F = a ji] = ka^/az is the Fresnel number. For large Fresnel numbers — )■ 




Figure 2. (color online) Analytical expressions and numerical evaluation for tr = 10 with a configuration 
average on 10 realizations as a function of atom number, (a) Results for sjq (blue circles), /jv (red squares) 
and analytical expressions 1/N + s^a (blue line), f^o (red line), (b) Forces acting on a cloud of atoms with 
Ao = —100 F: numerical evaluation (blue circles) of the average cooperative force. The full lines indicate 
(i) the force in presence of isotropic scattering, i.e. assuming /jv = (green dashed line), (ii) the force 
for continuous density distributions without disorder (red dotted line) and (iii) the total force taking into 
account cooperative scattering and disorder (blue line). 
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As illustrated by fig|2|3, the single-atom force is recovered in the limit of 
Nsoo ~ N/Aa'^ <C 1 i.e. for small optical thickness 60 ^ 3N/a'^. On the contrary, 
for large 69 the microscopic inhomogeneities can be neglected and cooperativity 
strongly modifies the radiation force. In particular, for small volumes the emis- 
sion is isotropic and foo ~ 0, whereas for large volumes the recoil at the emission 
compensates the recoil at absorption, /oo ~ Soo, which results in mainly forward 
emission. 



3. Experimental Results 

The analytical results presented in the previous section should apply to a large 
variety of clouds of cold atoms, including thermal cold atoms as well as degenerate 
quantum gases (as long as atom-atom interactions can be neglected). In this section 
we report first observations of cooperative scattering in different regimes : thermal 
clouds in a dipole trap or released from a dark MOT and for a Bose-Einstein 
condensate realized from a magnetic trap. 

3.1. Thermal Cold Atoms in a Dipole Trap 

A first series of experiments have been performed in Nice in 2008, where we were 
looking for possible signatures related to collective atomic recoil lasing (CARL), 
with spontaneous atomic bunching, when a large cloud of atoms is exposed to off 
resonant detuned light. In contrast to previous experiments (f71|) . this setup does 
not use a high finesse cavity, but the larger atom number we are able to trap 
() jgj) might allow to compensate for the absence of the cavity. The results of these 
studies did not shown any evidence of CARL, but provided the first signatures of 
cooperative scattering. The experiment has been performed using a vapor loaded 
magneto-optical trap (MOT) of *^Rb atoms. The cooling laser has been derived 
from an DFB master laser amplified by a tapered amplifier (TA), whereas the 
repumper laser has been a simple DFB laser. After a dark MOT period of 35ms, 
we load the atoms into a red detuned single beam dipole trap, formed by another 
DFB laser amplified by a TA and focused to a beam waist of ~ 200fim. The action 
of this dipole trap laser is twofold: on one side it holds atoms against gravity due 
to the dipole forces and on the other side, as the atom-laser detuning is not very 
large, the residual radiation pushes the atoms along the dipole trap. As one can 
see in Fig. [3] most of the atoms of the dark MOT are not loaded into the dipole 
trap and fall under the action of gravity. A small fraction of the atoms are however 
trapped in the dipole laser beam and are pushed along the axes of propagation 
of the laser. This observation is not surprising, given the moderate detuning (at 
least for dipole traps) we have used, with values ranging from Aq = —50GHz to 
Aq = —200GHz (detuning Aq given in respect to the F = 3 ^ F' = A transition of 
the D2 line of ^^Rb). This regime of dipole traps where radiation pressure cannot 
be neglected is not commonly studied, as spontaneous scattering of photons is 
usually not desired. The experimental protocol using the dipole trap to both hold 
the atoms against gravity and push them with off resonant radiation pressure made 
it difficult to independently change the parameters for the dipole trap (and thus 
the size and shape of the atomic cloud) and for the radiation pressure effects. This 
experiment however provided our first signatures of cooperative scattering which 
have been studied later on in a more quantitative way. We thus show in this section 
our first qualitative results, which could not be compared in a quantitative way 
to a theoretical model. As this experimental protocol required a finite interaction 
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time in the dipole trap in order to allow separation of the untrapped and trapped 
atoms, the shape of the clould of atoms changes during this interaction time. This 
further complicates a reliable comparison to theoretical model. 



1x10 2x10 3x10 4x10 5x10 

N . 



Figure 4. (color online) Average radiation force on the cloud of the atoms. For increasing atom number, 
the cloud is less displaced by the dipole laser. The interaction time for this experiment is 50ms, the laser 
beam power is WOmW for Aq = -76GHz (black square), Aq = -87GHz (red circles), Aq = -90GHz 
(green triangles) and Aq = —108GHz (blue stars). 



In Fig. [4] we plot the average force acting on the center of mass of the cloud of 
atoms kept in the dipole trap. This force can be extracted from the spatial dis- 
placement, as the interaction time is well known. The normalization of the force to 
the single atom force is roughly estimated from the measured values of laser power 
and detuning. The most striking point to notice in Fig. |4] is the clear reduction 
of the average radiation force with increasing atom number. This effect has been 
subsequently been studied in a quantitative way d^j). 

During the course of these experiments, we have observed some intriguing fea- 
tures which have not yet been studied in a quantitative way. For instance we found 
a systematic oscillation of the average radiation pressure force, as illustrated in 
Fig. [5j This feature clearly merits further experimental investigation, as if it is not 
due to an experimental artefact, it might be related to effects beyond the Markhov 
approximation, where ringing of superradiant time decay is predicted (0)- 
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Figure 5. (color online) Displacement of the center of mass of the atomic cloud as a function of detuning 
for different interaction times: 50ms (red stars), 60ms (open blue circles), 70ms (black squares) obtained 
for a laser power of P = lOOrnVF. 



3.2. Thermal Cold Atoms released from a dark MOT 

Following the first series of experiments and the development of our theoretical 
model presented above, we performed experiments which allowed for a quantitative 
comparison. The main results of these experiments have been presented in Q and 
we thus only show the main result here. These experiments have been done using 
the same vapor cell as in the previous experiments. These new results have however 
used the ^^Rb isotope. The experimental protocol used to allow for a quantitative 
measurement of cooperative scattering did not use atoms in a dipole trap, but 
atoms released from a dark MOT. Adapting the number of atoms interacting with 
the pushing laser by controlled repumping from the -F = 1 to the F = 2 hyperfine 
level in the ground state, we changed the optical thickness of the cloud without 
changing the size of the cloud. This allowed for a quantitative measurement of 
the radiation force as illustrated in Fig. [6j We take the good agreement between 
experiments and theory as a proof of the relevance of cooperative scattering in 
clouds of cold atoms. 




50 100 

bo 

Figure 6. (color online) Experimental data and fits using the cooperative radiation force (normalized to 
the single atom radiation force -Fi) in presence of disorder for Aq = — 1.9r (red squares) and Aq = — 4.2r 
(blue circles). The shadowed area corresponds to the non- physical region feo < 0. 
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3.3. Measurements in a magnetic trap 

The experiment described above present the first clean signatures of cooperative 
scattering of single photons along the lines of the model outlined in Refs. (0 [^l- 
At the large pump laser detunings and the large cloud volumes used, the light was 
strongly scattered into forward direction. In the case of small volumes, the absorbed 
light can be reemitted in directions other than the forward one. Cooperativity then 
also leaves its imprint in the geometry of the radiation pattern. 

Particularly small and dense samples can be made by cooling the atomic cloud to 
ultralow temperatures. Ultimately, their size is limited by the repulsive interatomic 
interaction making it difficult to reach sample sizes on the order of cr ~ 1. Another 
advantage of samples with temperatures below the recoil limit is, that the inter- 
action with a light field leaves detectable traces in their momentum distribution 
even when on average every atom scatters much less than a single photon. 

In an experiment performed at the university of Tiibingen we prepare a ®''Rb 
cloud in a magnetic trap and cool it down by forced evaporation to quantum 
degeneracy. The trap is ellipsoidal with frequencies coz/'^'^ = 25 Hz and cOrf^TT = 
160 Hz. We now apply in axial direction a light pulse with a power of Pq = 130 //W 
and a diameter of wq = 295 //m for a duration of tq = 10. ..500 /is detuned from 
the D2 line by Aq/ = ±0.5... ±4 GHz. Immediately after the pulse (within 100 /is) 
the trap is switched off. The cloud falls in free expansion for ttof = 20 ms before 
we apply an imaging pulse in transversal (radial) direction (see Fig. [7]). 




^^^P probe beam 



Figure 7. (color online) Scheme of the experiment. 

We either use thermal clouds of = 10^... 4 • 10® atoms at a temperature of 
T ~ 1 fiK or Bose-Einstein condensates of A^ = 10^. ..6 • 10^ atoms. In the case 
of thermal clouds the aspect ratio is 77 = 5.6 and the size is cr = 90, independent 
on A^ (see Fig. [s]). In the case of a condensate the interatomic interaction gives 
rise to a chemical potential of up to = 10 kHz and a transverse radius in the 
Thomas-Fermi limit of up to o" = 60. Also the aspect ratio depends on the atom 
number, because the mean field presses the condensate into the weakly confining 
dimension. Optical densities of up to 60 = 2000 are reached in axial direction. 

The effective radiation pressure is extracted from time-of- flight absorption images 
such as the one shown in Fig. [9} The first moment of the momentum distribution 
Ap = mAzcm/itof is a measure for the collective radiation force. The tendency of 
the radiation force to decrease with increasing atom number is clearly visible and 
will be subject to future investigations. 



4. Conclusion 



Cooperative effects in scattering of light by large clouds of cold atoms present phe- 
nomena which can be described by using driven timed Dicke states. These states 
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Figure 8. (color online) Blue lines and symbols correspond to T = 1 fiK cold thermal clouds (size and 
shape independent on A^), red ones to condensates in the Thomas-Fermi limit (size and shape depend on 
N). (a) Calculated (lines) and measured (symbols) A'-dependence of the radiation pressure ratio, (b) Same 
as (a), but in linear scales and as a function of resonant optical density. The experimental parameters were 
Ao = 500 MHz, (T+-pol, TO = 20 ^s, Tt^f = 20 ms. No parameter were adjusted. Every data point is an 
average over several measurements. 




Figure 9. (color online) Absorption images of atomic clouds released from the magnetic trap. 

present a convenient quantum approach even though the features exploited in the 
present paper do not go beyond what can be expected from a classical treatment. 
We have given detailed results on our theoretical model, with an analytical expres- 
sion for the modified radiation force when large clouds of atoms are used. Exper- 
imental confirmation of the modification of the radiation force has been observed 
in two different laboratories using different experimental configurations, ranging 
from cold atoms held in an optical dipole trap, atoms released from a dark MOT 
to atoms in a BEC setup with a magnetic trap. These experiments illustrate the 
wide range of situations where such cooperative scattering processes need to be 
considered. Important future extensions of this work arising from this coopera- 
tive processes would include any possible 'quantum' feature which could not be 
described with classical models. Fluctuations and higher order correlations seem 
appropriate first signatures to study in this respect. 
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